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Any theory that have a global spontaneously broken symmetry will imply the existence of very 
light neutral bosons or massless bosons (sometimes called Majorons). For most of these models we 
have neutrino-Majoron couplings, that appear as additional branching ratios in decays of mesons 
and leptons. Here we present an updated limits on the couplings between the electron, muon and 
tau neutrinos and Majorons. For such we analyze the possible effects of Majoron emission in both 
meson and lepton decays. In the latter we also include an analysis of the muon decay spectrum. 
Our results are \g £a \ 2 < 5.5 x 10 -6 , \g/j.a \ 2 < 4.5 x 10 -5 and \g r a\ 2 < 5.5 x 1CU 2 at 90 % C. L., 
where a = e, p, r. 

PACS numbers: 13.20.-v, 14.80.Mz, 13.35.-r 


I. INTRODUCTION 


Recently neutrino physics has given us many surprises with strong evidences for flavor neutrino conversion to 
another type of neutrinos. Analysis of data from solar, atmospheric and reactor neutrinos have shown us that no 
other mechanism can explain all the data, unless you have massive neutrinos jl, 2]. These experiments are the 
first strong evidence for non-conservation of family lepton number and this may indicate that new symmetries and 
interactions are the source of this phenomena. 

Experimental evidences of massive neutrinos imply that the Minimal Standard Model (SM) is no longer correct. 
The simplest extension would be the inclusion of right-handed sterile neutrinos, what would allow Dirac mass terms 
for neutrinos. Despite its simplicity, this approach does not help us to understand the neutrino mass scale or predict 
the neutrino masses. Due to the large gap between neutrino mass scales and the other SM scales, several mechanisms 
have been suggested to generate neutrino masses, relating this mass scale to new physics. In many of these models 
the masses are of the Majorana type or a mix between Majorana and Dirac types, what implies non-conservation of 
lepton number. 

As it’s well known, lepton number is an accidental global symmetry ({7^(1)) of the Standard Model. So, if the 
neutrino mass matrix includes Majorana terms, lepton number is broken either explicitly or spontaneously. If lepton 
number (L) is indeed a global symmetry 1 , its spontaneous breaking will generate a Goldstone boson, usually called 
Majoron Sjj, |4|. In this case the breaking of L sets a new scale and requires a scalar which carries lepton number and 
acquires a non-null vacuum expectation value (vev). Several extensions of the SM allow spontaneous L breaking and 
predicts the existence of the Majoron. However, the simplest extensions (with a triplet scalar) are excluded due to 
the experimental results of LEP on the Z° invisible decay. 


* Electronic address: andlessa@ifi.unicamp.br 
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1 In some grand unified theories (GUTs), lepton number is gauged and becomes a subgroup of a larger gauge symmetry. 
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Another important class of models which predicts the existence of Majorons are supersymmetric extensions of the 
Standard Model with spontaneous R parity breaking. In these models, the introduction of anti-neutrino superfields 
( N c ) and new singlet superfields ($) (which contain neutral leptonic scalars), allows spontaneous breaking of lepton 
number 5, 6, 7|. In almost all of these models the Majoron will be the imaginary part of some linear combination 
of sneutrinos, the scalar component of the super Higgs fields ( H u and Hd) and the $ superfields. Therefore we may 
safely assume as a mo del-independent coupling the following interaction term between J and v 2 : 


£ = ^2 ( 1 ) 

a,/3=e,/^,r 

where g a p is a general complex coupling matrix in the flavor basis. Because in most models J is basically a singlet 
(avoiding the constraints imposed by the LEP results), the above couplings are usually the most relevant ones to 
phenomenological analysis (at least at low energies). In most models we must also include couplings between neutrinos 
and a new light scalar (that we call y) with the same couplings as J: 


£ = J! 9a0V a V0X, ( 2 ) 

a,/3=e,/i,r 

Usually neutrino masses and mixings will depend on the vevs associated to the spontaneous breaking of L and the 
matrix g. In this context, knowledge of the couplings between neutrinos and Majorons may help us to understand 
the neutrino mass scale. However, this relation is very model dependent and may be very hard to realize in practice. 
Trying to make our results as model independent as possible, we will make no assumptions on g a p and present our 
results with and without the existence of the massive scalar y. Nevertheless, assuming Majorana neutrinos (what is 
reasonable since lepton number is violated), bounds on g a p may be transformed to the mass basis through the relation 

G = U T gU (3) 


where Gij is the neutrino-Majoron coupling matrix in the basis where the neutrino mass matrix is diagonal (M = 
diag(mi, m 2 , m 3 )) and U rotates the mass eigenstates to the flavor eigentates (see section UTTTI 1 . 

Majoron models can be interesting from the cosmological point of view because they can affect bounds to neutrino 
masses from large scale structure [8|. Neutrinos coming from astrophysical sources can also be significantly affected 
by fast decays, where the only mechanism not yet eliminated is due to neutrinos coupling to Majorons. This can 
affect the very high energy region, strongly changing the flavor ratios between different neutrino species f| or the 
lower energy region, as supernova neutrinos 0 , 11 ]- 

Presently we know that the role of neutrino-Majoron couplings is marginal in solar and atmospheric neutrinos, 
therefore it’s possible to put a limit on 12] 

'7 x 1CT 5 eU 2 


|G 21 | 2 < \g af 3 U* 2 U 01 \ 2 <4x10 


-6 


A 


( 4 ) 


2 The same being valid for non-supersymmetric models as well. 
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Category 

Upper Bound 

Process 

Reference 

solar neutrino constrain 

, r ,2 .pqn-e^xlO-'e^ 

|G211 <4Xl ° Am| ) 

V2 —*• J + v\ 

m 

supernova bounds 

\g ee \ < (1 - 20) x 10" 5 , 2 x 10” 11 < IdeJdMxl < 3 x KT 10 

v —> J + is, is + is —> J 

[10] 


lda/3 < 3 x 10 -7 or \g a0 \ > 2 x 10 -5 

ts—*J + is',is + is-^J 

[yj 

f3/30v decay 

|See| < 2 X lO" 4 

(A, Z ) —> (x4, Z + 2) + 2e + J 

m 

microwave background data 

Gij < 0.6110 -11 m^o and Gu < 10“ 7 

V —» J + v' 

M 

meson decay 

Ei=e,M,r l5ez| 2 < 3 X IQ- 5 , Ez=e, M ,r M* < 2 X 10" 4 

■k/K —> e + v + J 

[15, i6, W] 


TABLE I: Some of the previous bounds on neutrino-Majoron couplings from different sources. In the last two columns are 
shown the process used to constraint the couplings and the respective references. 


where G 21 is the neutrino-Majoron coupling in the mass basis and Atoq is the solar mass difference squared 
(Amlj = — m \). The observation of 1987A explosion ensures us that a large part of binding energy of su¬ 

pernova is released into neutrinos, what can be translated into the bounds [lol] 


Ideal < (1 - 20) x 10“ 5 , 2 x 10- 11 < \g eii \\g^\ < 3 x 10 


(5) 


such bounds where read off from Fig.l and Figs. 3,4 of Ref. m for Gee and \g e ii\\g 00 \, respectively. Also, limits from 


. Q 


decay and scattering of Majorons inside supernova give the bounds 


11 ] 


\g a p\ < 3 x 10 7 or \g a p\ > 2 x 10 % 


( 6 ) 


the first limit appears because if neutrino-Majoron coupling is strong enough the supernova energy is drained due 
Majoron emission and no explosion occurs; the second limit appears because if neutrino-Majoron coupling is too 
strong, the Majoron becomes trapped inside the supernova and no constraint is possible. 

While neutrinoless double beta decays (/3/30is) provide us the constraint 

|dee| < 2 x 10 -4 and \g ee \ < 1.5 (7) 

where the first (second) bound corresponds to Majorons with lepton number equal to L=0 (L=2) at 90 % C. L. 
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Also, no evidence of Majoron production was seen in pion and kaon decays and therefore QQuzi 

Y \9ei\ 2 < 3 x 1CT 5 , Y I^| 2 <2x10- 4 . (8) 

Besides the bounds mentioned above, there are bounds that depend on the rate of neutrino decay (v —> v' J). Such 
reaction depends on the neutrino lifetime, r, that is a function of neutrino-Majoron couplings in the mass basis, which 
we denote by G. Without additional assumptions on neutrino hierarchy, we can not relate directly the neutrino- 
Majoron couplings and the neutrino lifetime. One example is Ref. 14] that using cosmic microwave background data, 
put a stronger constrain 


Gin < 0.6110 


-11 —2 
^50 


and Gu <10 


-7 


(9) 


where 77150 = m/50 meV and G is the neutrino-Majoron in the mass basis: Gu and Gjj are respectively the diagonal 
and off-diagonal elements of G. 

Future experiments can improve the present bounds on many order of magnitude, we refer to Ref. 18, 191 for 
details. 

A summary of some of the previous bounds are shown in Table HI where we also show the respective relevant process 
used to constraint the neutrino-Majoron couplings. Almost all the bounds shown in TableJU assume one particular 
model or class of models. Probably the most model-independent result is from [l5j, [l6|, |17j, but in this case they 
assume not only neutrino-Majoron couplings but also neutrino-x couplings to compute the upper bounds shown in 
Table ID 

Here we will try to improve or make these limits more model-independent through an analysis of both meson and 
lepton decays. In Section III Al we discuss the limits from pion, kaons, D, D s and B decays, including decays of mesons 
into taus; in Sections III Bl and III Cl we include bounds from lepton decays (from the total rate and from the spectral 
distortions). We conclude transforming our bounds to the mass basis in Section Til Dl 

II. RESULTS 


Here we try to improve the bounds on neutrino-Majoron couplings through the analysis of possible effects on mesons 
and leptons decays as well as on the spectrum of the muon decay. We also rewrite our results in the mass basis, which 
in many cases is more important for model analysis. All the bounds obtained here have 90% C.L. and were obtained 
through the chi-square method assuming gaussian distributions and including both statistical and theoretical errors 
as follows 


2 (Rd ata - Rth eor ) 

<r ? . + a?, 

data theor 


( 10 ) 


where Rdata , Rtheor, o' data and o theor are respectively the experimental data of the rate R, the theoretical prediction 
for process R, assuming an incoherent sum of SM rate and Majoron contribution, the experimental error and the 
theoretical error. 
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A. Meson decay rates 


The process M —> l + vi was extensively studied in the literature and has the following total decay rate 20|: 


Tsm = 


Gl\V nn r\ 


87 r 


-frn m l m M 1 - 


M 


fradi 


( 11 ) 


where the f ra d accounts for radiative corrections. In Eq. CD. tom and mi are the meson and lepton masses, Gf is the 
Fermi constant, f m is the meson decay constant and V qq > is the respective Cabibbo-Kobayashi-Maskawa (CKM) matrix 


element. Unless specified otherwise we are using the quantities as listed in the Particle Data group compilation 


201 . 


We also use the same source to compute the relevant radiative corrections for the mesons decay rates. An important 
feature of Eq. m is that, because it’s a 2-body decay, Tsm is proportional to rrif , as it should be to conserve angular 
momentum. 

In the last few years several of the meson decay constants were calculated through lattice QCD 
used to obtain stronger theoretical predictions. We used both the experimental 


2 l|, which can be 


20 j and theoretical values 


22 , 


3, 24 


25j of f m on our calculations, but in most cases the results differ only by 10%. For this reason we will only show the 
results using the experimental values of f m . 

Due to the neutrino-Majoron couplings, the following process also contributes to mesons decay rates: 


M —> l + vi' + <7, 


( 12 ) 


where J stands for Majoron and may be any neutrino flavor. A complex analytic expression for the total decay 
rate is given in l.J . Here we show a simpler result valid in the limit mi = m„ = 0: 

G 2 F \V qq ' | 2 f 2 

m=e,fLT 


r, = 


768tt 3 


E \ 9lr] 


(13) 


This result shows that when Majorons are included, the total decay rate is no longer proportional to the lepton 
mass (since now we have a 3-body decay). Therefore, the Majoron contribution (Tj) may easily overcome the SM 
predictions (Tsm) if g ~ 1: 

Tj 


_ 1 —<L »1 

r SM 487 t 2 mf 

where we have assumed mi <C tttm- Assuming that the total decay rate is 


(14) 


T total — r SM + Tj, 


(15) 


the decay on J will be the dominant channel, unless g is small. Because only small deviations from the SM are allowed 
by experimental data, we must have g C 1. Following Eq. mi, we calculated upper bounds for \g a p \ 2 at 90 % C. L. 
The Table nn shows the bounds obtained through this analysis. As expected from the above remarks and the results 
on Table HT1 the most constrained matrix elements g a 0 will be those concerning e, since the approximation mi <C tom 
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Mesons 

\gea | 

ISea) 

1 g-ra\ 

Refs (exp. values) 

7T 

1.6 x 10” 4 

2.1 x 10" 1 


m 

I< 

9.5 x 1(T 4 

9.3 


m 

D 

1.6 x 10” 1 

2.3 

23 

[20] 

D s 


1 

6.3 

[26] 

B 

0.85 

1.5 

19 

[27] 


TABLE II: Upper bounds on X^= e a r 1 9ip\ 2 from meson decays with 90% C.L. The references for the experimental values used 
are shown in the last column. We only include the Majoron contribution, and not the new light scalar X- 


(16) 


is good in this case. We found that this bound can be improved using recent data 28] of the following ratio: 

rir+ ^ £+ +V e \ = (2 ' 416 ± °' 043) x 10 ” 5 
T(K+ —> + i/ M ) 

where the error is the quadrature of statistical and systematic errors. Because the Majoron contributions must be 
suppressed (as shown in Table HU, we may approximate the above ratio: 

r( /T+ -> e + + Ue ) T e SM + T e j ~ T e gM + T e j 

T(K+->n + + v,) r(( M + rrg M 5 1 j 

where T 6 ^ represents the decay rate with an e (/z) in the final state. In this way, using the previous statistical 
analysis, we can constraint the elements g ea (at 90% C.L.): 


Y \9el\ 2 < 1-1 X 10 

Z=e,/x,r 


-5 


(18) 


When it comes to the /i matrix elements (g nry ), the constraints from Table HT1 may also be improved if we consider the 
decay channels of mesons in four leptons [20]: 


BR(K + —> + v + u) < 6 x 10 


-6 


(19) 


Since the SM contribution to this decay is negligible, we may assume: 


BR(K+ +v v + J)< 6 x 10 


-6 


( 20 ) 


resulting on (at 90% C.L.): 


Y, M 2 <9xio- 5 


t=e,/i,r 


( 21 ) 

Finally, new experimental data for leptonic decay rates of heavy mesons such as the D + , D + s and B + mesons q,h 
271 . |29t l30| ]. allow us to impose limits to the r matrix elements ( g Ta ), as shown in Table HD The best bound being 


from the D+ leptonic decay on r + + v T (at 90% C.L.) 


26j|: 


Y |ffri| 2 <6.3 

Z=e,/x,r 


( 22 ) 


Because of large experimental uncertainty, this bound is quite weak, as can be seen above. 
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flQEzi 


We stress that unlike [15|, [lg, [13 the results shown so far do not include possible decays on a light scalar x and 
therefore are less model-dependent. If this new scalar is considered with a mass of 1 KeV (other choices for the x 
mass do not change these results as long as it is well below the initial state masses), the previous results are basically 
improved by a factor of 2 (again, at 90% C.L.): 


Y \9ea\ 2 < 5.5 x 10“ 6 , Y \ 9fla \ 2 < 4.5 x 10- 5 and Y lff™| 2 < 32 


(23) 


B. Lepton decay rates 


Because of its good experimental precision, lepton decays are good candidates for imposing bounds on neutrino- 
Majoron couplings. Moreover, in this case there aren’t uncertainties such as mesons decay constants and CKM 
elements. However, the leading term in T(li —> lj + 9j + Vi) is no longer proportional to the final lepton mass (as 
it was in the case of mesons), because the SM decay is a 3-body decay already. For this reason Tj < Tsm even for 
g ~ 1. In fact the inclusion of Majorons in the final state decreases the decay rate by a factor of s=s 10 (Tj ss Tsm/ 10, 
for g=l), instead of increasing it as it was in the case of mesons. Therefore we expect much weaker bounds in this 
case. But, as we will show below, it’s still possible to obtain good bounds for certain decays 3 . To calculate the 4-body 


[u, 32|. 


decay rate (T(Z —► l' + 9 + v + J)) we used the programs FeynArts and FormCalc 

As we did in the meson case, to constraint the g a p matrix we assume that the total lepton decay rate receives 
contributions from Majoron emission: 


r total(lct —■* 1/3 + + v a) — rj(/ Q —*l()-\-V + V + J)+ T SM^a ~^* 1/3 + v 0 + u a) (24) 

Because Majoron emission may change neutrino flavor (since g a /3 may be non-diagonal), Tj may have any type of 
neutrinos in its final state. For this reason we omitted the subindex in Fj. Besides, both neutrinos (v or D) may emit 
Majorons, what implies: 

r j{l a —*lp + i' + l'+J)oC ^^(Iffail 2 + \g/3s\ 2 ) (25) 

s 

where g a s and gps are the couplings between Majoron and the a anti-neutrino and /3 neutrino, respectively. In Eci. (I 25L 
the interference terms g a sg/3S are proportional to neutrino masses squared and were neglected. Because Table UTl shows 
that lighter leptons have stronger upper bounds, we will assume g a s gps- Therefore we will consider that Majoron 
emission by D a is dominant: 

rj(Z a Ip + 9 + v + J) oc Ifl'aisjP (26) 

S 


3 We thanks J. F. Beacom for suggestion to use lepton decays to constrain neutrino-Majoron decays. 










Using the experimental values for the p and r decay rates J2i 


and the same kind of analysis used in the last section, 


the following bounds were obtained at 90% C.L.: 

J2\9^\ 2 <4x 10- 4 , ^|.g rct | 2 <10xl0- 2 , (27) 

ol ol 

where the first bound comes from p decay and the second from r decay, both at 90% C.L.. For the r decay the same 
constraint is obtained if one considers decays in e’s or p’s. If we include the contributions from \ emission (again 
with mass of 1 KeV and at 90% C.L.): 

l-^«| 2 < 2.7 x 10- 4 , ^| 5 tq | 2 <5.5x10- 2 . (28) 

ol ol 

C. Spectrum of lepton decay with Majorons 


Another method that can be used to improve the limits obtained above is the analysis of the electron spectrum in 
the muon decay, which can be modified by the inclusion of Majorons. The normalized spectrum for the SM case and 
the Majoron case only are shown in Figure |T] 




FIG. 1: At left, normalized electron spectra for muon decay in the SM (solid line) and with Majoron emission only (dashed 
line). At right, the experimental allowed region (between solid lines) for the electron spectrum and the total predicted spectrum 
(SM plus Majorons) for three values of g 2 = |gp a | 2 . 


Precision measurements of the electron spectrum (used to impose constraints on non V-A interactions) may constrain 
g if we consider the changes in the SM spectrum after including Majoron emission. The usual analysis parametrizes 


the electron spectrum using two parameters (p and if) 


20 (: 


d,T(x) G 2 F rrvf 


, s 3 -a: 2 [3( 1 -a;) + ^p(4x-3) + 3)7 
48-7T 3 Euidx x 


m e 


1 — x. 


dx 


(29) 


where x = 


E 


Err 


m 2 + m 2 
and E max = —^- 

2 nin 


. For the SM the predicted values are p = 0.75 and ?y = 0: 


dTsAi(x ) _ G 2 F nVf 3 2 3 


dx 


4 8 „3 Ij* 1 -*! 


(30) 
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Previous Bounds 

Revised Bounds 

E a \g ea \ 2 <3x ltr 5 

£ q \g ea \ 2 < 5.5 x IO’ 6 

E a |<w| 2 <2.4x io- 4 

Iflwl < 4.5 x 10 

none 

£ q \g Ta \ 2 < 5.5 x 10” 2 


TABLE III: Comparison between the strongest bounds (including the scalar x) obtained here and the previous bounds from 
the same processes. All bounds are at 90% C.L. and the previous bounds are from m, mn. 


The current experimental values are p = 0.7509 ± 0.001 and g = 0.001 ± 0.024 20]. 
When the total spectrum (SM plus Majoron) is considered, we have found 


dr tota i(x ) G 2 F m 


dx 


[0-0066M 2 - 0m\g\ 2 x + (- + 0.35| 5 | V - (1 + 0.25| 5 | V] 


(31) 


where \g\ 2 = £ Q \g^ a \ 2 




FIG. 2: At left (at right), normalized electron neutrino (muon neutrino ) spectra for muon decay in the SM is the solid curve 
and with Majoron emission only is the dashed curve. In both cases we assume a diagonal g a g. 


From the above expression and Figure Q] we see that the most sensitive region is at the end of the spectrum (large 
a:), which can be used to constrain g. Figure [T] also shows the allowed region by experimental data (region between 
solid gray lines) and the shape of the total spectrum (including the SM and Majoron contributions) with different 
values of l5/xa| 2 - 

Because the spectrum is more sensitive to changes in the cubic term (or the p parameter), we consider the Majoron 
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contributions to p: 

Ptotai= |(2-0.25|g| 2 ) 

Using the chi-square method at 90% C.L. we obtain: 


E W <8 x 10- 3 


(32) 


(33) 


As can be seen in Figure[2l the Majoron main modifications to the spectra occurs in the neutrino spectrum, which 
has been measured by the Karmen experiment { 33 1 . However, due to experimental uncertainties, the resulting bounds 
on g are too weak in this case. 

Summarizing, the strongest bounds are given in the Table mu where we compared the previous limits and the 
newest constraints obtained here. 

All bounds from Eqs. G51) . KS1) . (1351) can be written as 

E \9io\ 2 <Lt (34) 

o:=e,/z,r 

where Lj is the strongest upper bound for Y2 a \ 9 ia \ 2 (see Tabic Hill) . From this constraints, we assume the conservative 
limit, where the upper bound applies not only for the sum, |<?z«| 2 , but also for the individual elements, as | gi a \: 


\fjia | < L h V a = e,p,r 


(35) 


D. Mass Basis 


All the results obtained so far are written in the flavor basis. However, in many cases, theoretical analysis are easier 
on the mass basis. We have two possible cases: Dirac or Majorana neutrinos. In this section we assume Majorana 
neutrinos to transform our bounds to the the mass basis. 

We can translate the previous results to the mass basis using the transformation matrix U 


20 |: 


/ 


U = 


n —iS 


C12C13 S12C13 si 3 e 

—S12C23 _ ci2S23Si3e I<5 C12C23 — Si2S23Si3e l<5 S23C13 I x diag(e lai ^ 2 , e*“ 2 ^ 2 ,1) 


(36) 


\ S12S23 — Ci2C23Si3e I<5 —C12S23 — Si2C23Sl3e i<5 C23C13 


where Cy = cos(0y) and = sin( 6 ij). The neutrino mass matrix is given by M = diag(mi, m 2 , m 3 ) and for a given 
mass mi, we can written all other masses as a function of mi and the squared mass differences as follows 


A mf 2 = m-2 ~ m i = and = m3 — = A m 2 atrn . 


(37) 


Although the mass differences and angles have been measured experimentally 


20| |. we have no information on the 


Majorana phases <5, aq and aq. To calculate the bounds in the mass basis, we will use the transformation rule for 
Majorana neutrinos 


G = U T gU 


(38) 
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FIG. 3: Exclusion curves in the m i — G plane. We show the curves for gi, g 2 and g$. These curves are for #13 = 0, non-zero 
values of 6*13 are more restrictive. 

where g is the neutrino-Majoron coupling matrix in the flavor basis and G is the neutrino-Majoron coupling matrix 
in the mass basis. 

Although it is not valid in general, many models 


M 


35 


36] have the following property (at least in some limit) 

G = diag(g 1 ,g 2 ,g 3 ) oc M = diag(rn 1 ,m 2 ,m 3 ) (39) 


Following 


37| . we calculate the allowed region for different values of <5, a\ and a 2 and then choose the union of 


these regions as the final result, valid for any value of the phases, as shown in Figure [3) 


III. CONCLUSIONS 


Using three different techniques we were able to constraint the neutrino-Majoron couplings. The strongest con¬ 
straints are shown in Table HTT1 Considering only the limits from meson decays we improve by one order of magnitude 
the previous limits on |< 7 ect | 2 and \g^, a \ 2 }15l . 16, 17]. Although the best constraints were obtained from meson decay 
rates, we have shown that independent bounds can also be obtained from g and r decays. The latter one being 
the best to constraint the g Ta elements. We stress that the bounds on g Ta shown in Table mu is probably the first 
model-independent constraint for this parameter. 

The third alternative used was an analysis of the spectrum of muon decay. Despite its potential for constraining the 


S/xct elements, the experimental values are not precise enough to make such an analysis useful. Our best constraints 
are \g e a\ 2 < 5.5 x 10“ 6 , Ig^a] 2 < 4.5 x 10 -5 and \g T a\ 2 < 5.5 x 10 -2 , a = e, g, r at 90 % C. L. . 
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Because the models cited here usually try to explain the neutrino mass scale, it may be convenient to analyze the 
limits on neutrino-Majoron couplings in the mass basis. With that in mind we transformed all our results from the 
flavor basis to the mass basis, using the current values for the angles of the neutrino mixing matrix. As shown in 
Figure [3] the constraints on the mass basis are usually weaker than those on the flavor basis. 
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